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The Hamiltonian control of n qubits requires precision control of both the strength and timing of 
interactions. Compensation pulses relax the precision requirements by reducing unknown but sys- 
tematic errors. Using composite pulse techniques designed for single qubits, we show that systematic 
errors for n qubit systems can be corrected to arbitrary accuracy given either two non-commuting 
control Hamiltonians with identical systematic errors or one error-free control Hamiltonian. We also 
examine composite pulses in the context of quantum computers controlled by two-qubit interactions. 
For quantum computers based on the XY interaction, single-qubit composite pulse sequences nat- 
urally correct systematic errors. For quantum computers based on the Heisenberg or exchange 
interaction, the composite pulse sequences reduce the logical single-qubit gate errors but increase 
the errors for logical two-qubit gates. 

PACS numbers: 03.67.Pp; 32.80.Qk; 82.56.Jn 



I. INTRODUCTION 



The control of quantum bits for quantum computation requires a high degree of accuracy. Aside from coherence 
and random noise, systematic errors limit our ability to control these quantum systems. These errors include slow 
fluctuations in control parameters relative to the experimental time and slight imperfections in fabrication. Overcom- 
ing these systematic errors will be crucial to achieve the potentially high-accuracy gates required for fault-tolerant 
quantum computation with local gates [J, 0, H[ . The problem of unknown systematic errors has been studied exten- 
sively in NMR 4]. In NMR a large collection of spins are addressed by an RF field with an unknown spatial variation. 
To overcome this variation, broadband composites pulses were introduced [H, [fij]. 
, In principle, compensating pulses can be used to correct unknown systematic errors in single qubit gates to arbitrary 
■ order Q. In a real experimental situation, other errors begin to accumulate and higher-order pulses may be of limited 
use @. The second order broadband pulse devised by Wimperis (BB1) Q is the standard of compensation and has 
been extended to two-qubit couplings by Jones (Toj . In this paper, compensation pulses for multi-qubit systems and 
Hamiltonians are examined using BB1 as an example pulse. BB1 and the higher-order pulse sequences of @, HH are 
fully compensating; the pulses do not require a specific input state of the system and can be used to replace single 
0© pulses that are part of a larger sequence. 

The paper is organised as follows: Section [TT] describes how the control theory (and related geometry) of multiple 
qubits is suited for the type of compensation pulses used on single qubits. Section [Ml introduces our notation and a 
generalised BB1 sequence. Section llVI reexamines the two-qubit pulse sequence of Jones 10] in the case of multiple 
systematic errors. Different methods for creating BBl-style sequences are compared. Section IVl generalises to n qubits 
and proves inductively that only two systematic errors need to be correlated to achieve arbitrary correction in all 
systematic errors. Section IVII examines cases with sufficient control for universal quantum computation but not full 



control of the n qubit space. Finally, we conclude in Section fVIH 



II. CONTROL THEORY AND GEOMETRY OF n QUBITS 



The model we consider is n qubits and M dimensionless Hamiltonians denoted H m . We define a pulse as applying 
H m with constant strengths f2 m for a time t where fl m is bound between — f2 max and f2 max - The resulting unitary 
evolution is U{t) = exp(— i^2 m £l m H m t). The applied pulse may not create the desired evolution due to systematic 
errors in the control strength Cl' m — £l m (l+5 m ) and the timing t' = t(l+St). In this model, timing errors are correlated, 
while the individual strengths could have independent errors. The source of the errors will not be considered and we 
will examine unitaries of the form £/({#}, {e}) = exp(— i J2 m 8 m (l + e m )H m ). 
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The quantum system is universally controllable without unknown errors if H m generates the entire control algebra 
of su(2 n ) by addition and the Lie bracket [l2j. The very same technique can be used to determine if a composite pulse 
sequence exists [l3j]. Additionally, the Lie bracket can be used to constructively build pulses, e.g. the Solovay-Kitaev 
composite pulse sequences in Q- 

For n qubits the corresponding Lie Algebra is su(2 n ). We choose as a convenient representation of the generators 
of the algebra, rjj = i (SC=i mod 4 fe )/4 fc -!j where <tq = I is the identity on the qubit and a\ = X, cr 2 = Y, and 
03 = Z are the single qubit Pauli operators. 

There are 4™ — 1 operators since the generator of the global phase 7/0 = \ ®fe=i I i s outside of the algebra of su(2 ra ). 
For any two generators r\i and 77^ , we find that either they commute [rii,r)j]=Q or [r)i,r)j] — ieijkijk- If they do not 
commute, the two operators generate a representation of su(2). 

The Lie algebra then imposes that given Pauli-operator generators with the same systematic control error, arbitrarily 
accurate composite pulses can be created, if and only if they do not commute. Furthermore, if they do not commute 
the resulting pulse sequence will have the same form as a single qubit pulse sequence [131 ] . A geometrical interpretation 
is that controlling two elements that do not commute is homomorphic to rotations on a sphere while the space for 
commuting elements is a 2-torus [l4|, EEI • 



III. NOTATION AND BB1 REVISITED 



The goal is to create accurate multi-qubit unitaries in the presence of systematic errors in 9i . For each case, we will 
start by defining the set of generators we control, {Hi}, and denote the unitary transformations as 

Ui(9i) = 

Ui >m (6i,d m ) = exp(-i(6iHi + 6 m H m )) 
U^ mi n(0i,0 m ,9 n ) — exp (—i(6iHi + 6 m H m + 9 n H n )) 

Ul y m v n,p{dl,0 m ,9 n ,0p) = .... (1) 

We will be particularly interested in sets of three generators that have commutation relations equivalent to su(2). 
In this case, rotations around the sphere can be used to guide the mathematics. The compensation pulses we present 
require that we can perform both positive and negative rotations. Physically this corresponds to inverting applied 
fields and changing the sign of multi-qubit interactions. 

A useful metric for evaluating the effects of control errors is the infidelity, 1 — F(U, V), where F is the fidelity, 

F(U, V) = min J WV (tp\ VW \tp), (2) 

where U is the ideal unitary and V is the actual operation affected by the systematic error e. We choose this 
measurement over the distance, D(U, V) = \\U — V\\, to avoid complications due to a global phase, e.g., U — X and 
V = —X. For U — Ui(6 1 ) and V — + ei)), the distance scales as 0(e) and the infidelity scales as 0(e 2 ) 

EE!. 

Imagine we would like to perform U\{9) but our systematic control errors limit us to control of the form Lq ; 2(#i(l + 
#2(1 + £2))- Compensation sequences minimise the effect of these errors by applying successive error-prone pulses 
that cancel the leading error terms. In this notation, the BB1 sequence Q is 

V w {9, Hx, H 3 ) = U X {6{1 + ei))2V(& H U H 2 ), (3) 

where Tw(4>, H\,Hq) is the correction sequence with <j> = acos(— Q/Att) 

T w ((f > ,H 1 ,H 2 ) - Ui,2(7rcos(^)(l + e 1 ),7rsin(0)(l + e2)) 
xJ7i.2(2ttcos(30)(1 + ei), 2tt sin(3</>)(l + e 2 )) 

x[/i,2(ttcos(0)(1 + ei),7rsin(0)(l + e 2 )). (4) 

We refer to this sequence as BB1-W and when t\ = e 2 = e the sequence yields an infidelity that scales as e 6 , 
1 — F(Vw(0, Hi, H2), Ui(6)) = 0(e 6 ), or a distance that scales as e 3 , details in [A] An infidelity that scales as e 2n 
corresponds to a distance that scales as e™ [3] • The fine control of the relative amplitude or phase cf> allows for the 
correction; the compensation of higher order terms relics on increasingly finer control. 

The BB1 pulse sequence was derived in the context of single spins in NMR where Hi = ^X and Hi = |Y [9(. In 
many controlled quantum systems, the control occurs in a rotating frame and the difference between applying the 
generator Hi or cos((f))Hi + sin((f>)H2 is phase shifting the applied oscillating field relative to the rotating frame [HI]. 
As a result, for single qubit gates it is often reasonable to assume ei = e 2 . 
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IV. TWO QUBITS AND MULTIPLE ERRORS 

Jones applied BB1 to two qubit gates jioj]- His construction assumes that the single qubit gates are without error. 
In the context of NMR, the natural two-qubit Hamiltonian is Hi = ^ZiZ 2 . The error in the control of H\ is unrelated 
to the error in H 2 — \X\, in this case no error. The direct application of BB1-W by simultaneous Hi and H 2 pulses 
would fail to correct the errors. However, the error free rotations about Xi allows us to construct unitaries that are 
generated by H% — ^Y X Z 2 . Since the algebra is equivalent to rotations, we can use a y (H 2 ) rotation to rotate the x 
axis (Hi) to an axis in the x — z plane (ciHi + c 2 H 2 ) yielding 

U 2 {^)Ui{9)U 2 {-4>) = l7 li3 (0cof#),0siiO)). (5) 

This identity was used by Jones [l(| to create an alternative pulse sequence, we will refer to as BB1-J. BB1-J 
transforms the requirement of relative amplitude-control (BB1-W) into the accurate control of a rotation. We note 
that in NMR the sign of the ZZ Hamiltonian Hi is determined by the molecule [f| • This shows that not all of the 
control Hamiltonians require invertible couplings in order to compensate. The correction sequence is then 

Vj(6,Hi,H 2 ) = + ei))U 2 (cf>(l + e 2 ))*7i(7r(l + ei))tf 2 (4>(l + e a )) 

x*7 2 (30(l + e 2 ))U 1 (2n(l + e 1 ))i7 2 t (30(l + e 2 )) 

*U 2 (0(1 + e 2 ))Ui{n(l + ei))Ul(4>(l + e 2 )). (6) 

This sequence yields an infidelity that scales as ef when €2 = [10|]. The scaling for when e 2 ^ is examined inlAl 

The utility of the any fully-compensating pulse sequence is that it can be used to replace single pulses in a sequence. 
If Xi and Yi have the same systematic error, we can correct the Xi rotation by BB1-W before correcting the ZiZ 2 
transformation by BB1-J. The sequence of BB1-WJ is 

V WJ {9, H U H 2 ,H 4 ) = Ui{6{\ + ei))VW(0, H 2 , H 4 )Ui(tt(1 + et))V w (<P, H 2 ,H 4 ) 

x V w (30, H 2 , Hi ) Ui (2tt (1 + e 1 ) ) V w (30, H 2 , ff 4 ) 

x V w (cf>, H 2i H 4 )Ui(tt(1 + ei))V w (<f>, H 2 ,H 4 ), (7) 

where 

V W {4>,H X ,H 2 ) = Tl HiH2 ( a co S (~^/47r))ulm + ei)). (8) 

This sequence replaces error prone U 2 (<j>(l + e 2 )) pulse with the corrected rotation Vw generated by the BB1-W 
sequence. Here, {Hi, H 2 , H 3 = — i[Hi,H 2 ]} is a representation of su(2) and {H 2 ,H4 — ^Y,H 5 = — i[H 2 ,Hi\} is also 
a representation of su(2). The assumption is that the errors of H 2 and H4 are equivalent, e 2 — 64. The infidelity then 
scales as (aef + (3eie 2 ) 2 where a and j3 are constants that depend on 9, Hi, H 2 , and H 4 . 

For fixed e 2 , the infidelity at small ei scales as el in e\. This is the same order as the uncorrected pulse in ei, 
although with a substantially smaller infidelity. In the case of Hi = \Z X Z 2 , H 2 = \X t , and H 4 = \Y X , where 
e 2 = ex — 0.01, the infidelity in this regime is a factor of 10 8 smaller than the uncorrected pulse (see Figure [l}. 
For Vwj, the infidelity scales as ef when e x < ^iE 2 2 . However, we can replace the BB1 — W sequences Vw in Vwj 
with higher order pulse sequences, for example the Bn sequences where B2=BB1 Q. In this case, the infidelity will 

scale as (ae\ + 7„eie 2 " +1 ^) 2 , where 7„ is a constant that depends on 9 and Bn. As a result, the value of e x where 
the scaling changes from e\ to t\ becomes smaller and smaller. In Figure [TJ we compare the scaling properties of the 
BB1-WJ and the higher order BB1-WJ where we have replaced the Vw BB1 sequence with the B4 sequence ^ s 
expected, the error ei = tzz where the scaling changes from e 6 to e 2 changes from w 10~ 2 for BB1-WJ to w 10~ 4 for 
BB1-WJ. In principle, given a target infidelity and systematic errors e < 1 [7[, we can construct a pulse sequence with 
an infidelity guaranteed below the target infidelity. We note that in practice other errors including random control 
errors and decoherence typically limit the fidelity. 

These sequences are each optimised for different correlations in the errors. BB1-W performs well when errors in the 
control of ZiZ 2 and Xi are correlated while BB1-J is optimised for when one control has no error. BB1-WJ combines 
both strategies by first correcting the correlated errors and then correcting the independent error. 

In Figure^ we compare the ideal unitary U — Uzz(k /4) = exp(— i~ZiZ 2 ) to the approximate unitaries V assuming 
errors equivalent errors in X\,Y\ and uncorrelated errors in Z X Z 2 . BB1-J (V = Vj(n/4,ZiZ 2 /2,Xi/2)) outperforms 
BB1-W (V = VV( 7r /4, ZiZ 2 /2, Xi/2)) when either error is low. BB1-W is preferable when the systematic errors are 
identical. BB1-WJ (V = Vwj{k/4, ZiZ 2 /2, Xi/2, Yi/2)) results in low errors over the range of two errors. Initial 
compensation of the Xi pulses results in better compensation of ZiZ 2 . 
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FIG. 1: Comparison of BB1-WJ and the higher order BB1-WJ pulse sequences applied to a £/zz(7r/4) operation. For a fixed 
Xi and Y\ error ex, the infidelity after a BBf-WJ correction scales as (ae% z + f3ezze x ) 2 (see text). For the same ex the 
BB1-WJ sequence scales as {ae 3 zz + "f4ezze x ) 2 , extending the regime where the infidelity scales as e% z - 
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FIG. 2: Comparison of (a) BBI-W, (b) BBf-J, and (c) BBf-WJ pulse sequences applied on Uzz(tt/4:) operation on a pair of 
qubits. BBf-WJ assumes X\ and Y\ have equivalent systematic errors. 



V. EXTENSION TO MANY QUBITS 



Given a control operator with a systematic error and a perfect rotation that transforms that operator to an 
orthogonal independent operator, we can perform compensation, e.g. BB1-J. Given two control operators with 
correlated errors that are generators of su(2), we can perform compensation, e.g. BBI-W. As a result, in principle 
one can perform arbitrarily accurate composite pulses on a controllable quantum system where all the controls have 
independent errors except two. 

As an example, imagine n qubits in a row with single qubit operators and tunable Ising couplings. The Hamiltonians 
are Xj, Yj on each qubit and ZjZj+% between neighbours. If for the qubit n, X n and Y n have uncorrelated error, there 
does not exist a compensation pulse 13]. However, if the X and Y systematic errors are correlated on the the first 
qubit but otherwise independent, the following sequence can be used to generate an arbitrarily accurate X rotation 
on the nth qubit. 

For the initial qubit with correlated X\ and Y\ errors, BBI-W is used. To correct Z\Zi, BB1-J is used with 
BBI-W corrected X\ pulses. This is the sequence BB1-WJ. X 2 on the second qubit is then corrected via BB1-J using 
BB1-WJ corrected Z\Z 2 pulses. We denote this sequence as BB1-WJJ or BB1-WJ 2 . Errors on the nth qubit can be 
compensated by repeated use of BB1-J along the chain, first correcting Xj, then ZjZj + \ and then Xj+i until X n is 
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Error 



FIG. 3: Compensation of £7x„(7r/4) by application of BBl-WJ 2 '" -1 ' . Compensation of Ux n by BB1-W pulses using Y n works 
only when the errors are correlated. Anticorrelated errors between X n and Y n increase the infidelity. BBl-WJ 2 '™ -1 - 1 uses the 
correlated errors of X\ and Y\ and a chain of ZjZj+\ interactions to compensate the X n rotation. The results for X2, X4 and 
X§ are shown. 



reached. The total sequence correcting the nth X rotation is denoted BBl-WJ 2 '™" 1 -'. 

Figure [3] compares correcting a 7r/4 X rotation as a function of chain length assuming equal magnitude errors 
for all operators but with a random sign except for X\ and Y\. The correlated and anti-correlated lines serve as 
references. If X n and Y n have correlated errors, then local BB1-W greatly reduces the infidelity. In the worst case 
scenario, the errors are anticorrelated and the compensation pulses add additional error to the initial overrotation. 
X n rotations can still be corrected using BBl-WJ 2 *™ -1 ), if only X\ and Y\ are correlated. The error increases with 
position (comparing BB1-WJ 2 to BB1-WJ 10 ) on the chain for large errors but approaches an equivalent fidelity for 
small errors. Asymptotically, the correction of X n rotations by sequential correction (BBl-WJ 2 *" -1 )) is equivalent to 
the BB1-W correction composed of correlated X n and Y n rotations. Replacing BB1 with the pulse sequences from Q 
allows for the creation of arbitrarily accurate pulse sequences. 

Although, this is not practical on a large scale, it can lead to a constant reduction in the number of gates that need 
to be calibrated at the beginning of an experiment for a large quantum system. Per region of computation, only a 
few highly reliable quantum gates can be used to reduce systematic errors in their neighbours. 



VI. LIMITED UNIVERSALITY 



An interesting theoretical proposal with potential applications for quantum dots \19L l2fj , superconducting qubits 
PH , and trapped ions [22j is the use of only two-qubit interactions for quantum computation [23[ . These two-qubit 
interactions are chosen to generate a sufficiently large algebra to create universal computation on a subspace of the 
total Hilbert space. We examine composite pulses for XY and Heisenberg interaction based quantum computers. 
The pulse sequences require that the sign of the two-qubit interactions can be inverted. Coupled quantum dots have 
been shown to exhibit reversible exchange couplings which can be controlled by an external magnetic field [25l . |26| , 
and may be promising candidates for this encoding. 



A. XY 



A Hamiltonian made of XY interactions, -*4{ij} = ^(XiXj + YYj), has been shown to be universal over 3-qubits 
encoded into one-qubit [H,[24|]. The XY interaction preserves the projection of angular momentum along z but does 
not preserve total angular momentum. The qubit is encoded in a subspace of the qutrit defined by m z = 1/2 or 
m z = -1/2. 
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FIG. 4: Application of BB1-W pulse sequence correction on an effective ?7z(7r/4) operation (T gate) on the XY qubit following 
[24| . Z rotations are implemented with the five pulse V3 sequence. Each pulse in the sequence is replaced with a BB1-W 
compensating pulse sequence. 

Irrespective of the encoding, compensation is possible if the systematic errors are shared because 

A {1 , 2} = ^(X 1 X 2 +Y 1 Y 2 ) (9) 

^{2,3} = ~(X 2 X 3 +Y 2 Y 3 ) (10) 

A'= [A {h2} ,A {2j3} ]/(i) = ]^{X X Z 2 Y 3 -Y X Z 2 X 3 ) (11) 

is a representation of su(2). For three qubits, we can block diagonalise the operators into four irreducible repre- 
sentations. For m z = ±3/2, the irreducible representation is one-dimensional. For m z = ±1/2 the irreducible 
representation is three-dimensional. 

The operators X, Y and Z act as the Pauli matrices on the encoded space and can be performed with pulses 
utilising the XY interaction alone. For example, Z rotations on the encoded qubit are implemented by the five pulse 
V3 sequence which uses XY interactions between each of the physical qubit pairs 0, H3| • 

W, e) = C/{i,2}(-f (1 + e))tf {2 ,3}(-§ (1 + e))C/{i,3}(-f (1 + e)) 

xC/ {2i3} (f(l + e))C/ {1)2} (f(l + e)) (12) 

When e = 0, V 3 {9, e) is equivalent to a rotation about the z axis in the code space, Uz{0) 1 up to a global phase. The 
V3 sequence explicitly requires invertible couplings between physical qubits, which may limit the types of systems 
that an XY computer can be built from. Assuming the errors are proportional for each H^ijy, we can correct the 
timing error using BB1-W for each pulse; each Ui k jy(9(l + e)) is replaced by Vw{0, Ar k j^,A^j ^k})- The results of 
using the correction are shown in Figured] 

The remarkable part of the XY interaction is that the su(2) algebra of neighbouring XY operators is independent 
of our choice of encoded qubit. The exact same methods can be used to compensate the two-qubit gate sequences. 
Furthermore, we can apply our results for the ZZ chain from Section [V] to show that only two neighbouring XY 
interactions need to have identical systematic errors. 



B. Heisenberg 

The Heisenberg Hamiltonian G{ij} — (XiXj + YiYj + ZiZj) has also been shown to be universal [H, [29|, [30j ] . 
Furthermore, for certain arrangements of spins and exchanges it serves to protect errors by both energetics and 
symmetries [31(. It is more convenient to write this as j\ = {—I + 2E(i,j)) where E(i,j) exchanges the states of 
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FIG. 5: Application of BB1-W pulse sequence correction for an l/zi^/A) operation on a Heisenberg qubit encoded on three 
spins. The qubit is encoded on a su(2) subspace on which compensation pulses are performed. BB1 sequences greatly improve 
fidelity for states in the code space but fail for states with support outside the code space. 



qubits i and j. As a result [Guj\, Gy^}] — j, k) — P{i 1 k,j)) where P(i,j, k) is the cyclic permutation of i, j, 

and k. This does not result in a representation of su(2). 

The exchange Hamiltonian E(i, j) preserves both total angular momentum and the projection along z. For a single 
logical qubit made of three spins [30], any two exchange terms represent the algebra u(l)@su(2). The qubit is 
encoded into the su(2) block corresponding to the total angular momentum, S = 1/2, and projection, m z = 1/2. The 
exchange Hamiltonian generates rotations equivalent to su(2) on the code space with E(l, 2) and -^(E(l, 2)+2E(2, 3)) 

corresponding to Z and X respectively. Up to a global phase, G{ ij } and 2E(i,j) generate the same unitary evolution. 
Assuming that G{ 12 } and G{ 2 ,3} have equivalent systematic errors and the interaction strengths can change sign, 
compensation is then possible on the code space using BB1-W. 

It is not clear that the two-qubit gates can be corrected since the state has support on the u(l) and su(2) blocks. 
In Figure [5] we calculate the three-qubit fidelity after compensation. If we limit ourselves to states in the code space, 
the compensation works as expected. Allowing states outside of the code space, the compensating pulses are worse 
than the uncompensated pulse for low errors. This result is expected since the operators do not form an su(2) algebra 
which the BB1 sequence depends on. 



VII. CONCLUSIONS 



We have shown that arbitrarily accurate compensation is possible with a fully controllable system if either two 
non-commuting Hamiltonians that generate su(2) have equivalent systematic errors or if a single Hamiltonian is error 
free. In the case of two non-commuting Hamiltonians with equivalent systematic errors, pulses of both positive and 
negative amplitudes are required. The underlying pulse sequences are equivalent to sequences for qubits. Furthermore, 
we reemphasise the importance of the algebra and show that the same compensation pulses work for universal XY 
quantum computation but not for universal Heisenberg quantum computation with a three-qubit encoding. 

Compensation pulses are well-suited for single qubits controlled by interaction with electromagnetic waves in the 
rotating frame. In this case, the difference between X and Y Hamiltonians is simply a change in the phase of the 
electromagnetic wave. Uncertainty in the amplitude of the applied field naturally leads to an unknown but equivalent 
error in the two Hamiltonians. 

For systems based on two-qubit interactions, the requirement of positive and negative couplings can present a 
challenge. Furthermore for solid-state systems, the gate couplings will most likely have independent systematic 
errors. Our work presents a possible solution. For the XY case a single, invertible error-free XY interaction can be 
used to build accurate gates. Using the identity that exp(— iirA{i 2 })A{2,3} exp(i7rA{ 12 }) = — A{ 2 .3}, we can create 
effective negative couplings for neighbouring XY interactions. We can then generate arbitrarily accurate unitary gates 
locally using BB1-J type sequences. These can interact with their neighbours, etc. This is impractical but it does 
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suggest that a system with a few low-error invertible couplings could efficiently compensate neighbouring high-error 
couplings. 

The su(2) algebra underlying these compensating pulse provides additional incentive to continue development of 
single qubit compensation pulses. Shaped pulse sequences or continuous time control can lead to further improvements 
[32J. The question remains how to develop composite pulses that do not rely on a su(2) or so(3) subalgebra. The 
Lie algebraic technique of 13] rules out composite pulses with Heisenberg coupling. However, we know that over 
an encoded space at least space single qubit compensation pulses are possible. The development of compensation 
pulses that do not use the geometry of the sphere and the development of techniques for identifying compensation 
compatible subspaces are both interesting challenges. 
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APPENDIX A: ANALYTICAL EVALUATION OF ERRORS IN BB1 AND BB1-J 

The original BBf sequence @ is BBf-W where both errors are equivalent 

V{9, H U H 2 ) = E/i(0(l + e))T{<j>, H U H 3 ), (Al) 

where T(cf>, Hi, H 2 ) is the correction sequence with <f> = acos(— 9/4tt) 

T(4>, Hi,H 2 ) = C/ 1 , 2 ( 7 rcos(0)(f + e),7rsin(^)(l + e)) 
xf7i )2 (27rcos(3^)(l + e),27rsin(30)(l + e)) 

x U 1>2 (tt cos(0) (1 + e) , 7T sin(0)(l + e)) . (A2) 

The 7r rotations toggle 3<p to -<j) (see Q). After removing identities, we can rewrite T((f>, Hi, H 2 ) as 

T{(j),Hi,H 2 ) = C/i, 2 (7recos(0),7resin( ( /.)) 
xt r i, 2 (27re cos(0), — 27re sin(</>)) 

X 1/1,2 (7recos(^), 7resin(0)). (A3) 

We use the Magnus expansion 33] to combine the three unitary operators, 

T{cj),Hi,H 2 ) = cxp(ie6Hi+ie 3 M 3 (cb,Hi,H 2 ) + 0(e 5 )) 

= Ui(~e0)[l + ie 3 M 3 {0,Hi,H 2 ) + O(e 4 )] (A4) 

where 

2tt 3 

M 3 {<j), Hi,H 2 ) = — cos(<?!)) sin 2 (0)i/i + 2tt 3 cos 2 (</>) sm{(j))H 2 . (A5) 

The second order term vanishes due to the symmetry of the pulse sequence 0,0- As a result, 

V{9,Hi,H 2 ) = [7i(0(l + £ ))U 1 (-e6) [l + ieHd 3 (4>,Hi,H 2 ) + 0{e A )] 

= Ui{6) [l + ie 3 M 3 (0,Hi,H 2 ) + 0(e 4 )] . (A6) 

This shows that V(9,Hi,H 2 ) is an approximation of Ui{6) that scales as e 3 in distance [?J and e 6 in infidelity [To! ]. 
Both the distance and infidelity depend on the specific Hi and H 2 We perform a similar analysis for BBf-J. Starting 
from Equation [6l we find 

Vj(9,Hi,H 2 ) = Ui(8(l + ei))Tj(<S>,Hi,H 2 ), (A7) 

where Tj($, Hi, H 2 ) is the correction sequence with $ = acos(6*/47r)(l + e 2 ) = cf>(l + e 2 ) 

Tj{^,Hi,H 2 ) = C/i, 3 (7rcos($)(f + ei),-7rsin($)(l + ei )) 
xC/ 1 , 3 (27TCOs(3$)(l + ei), -27rsin(3$)(l + ei)) 

x[/i, 3 (7rcos($)(l + ei), -7rsin($)(l + e)i), (A8) 
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with H$ = i{Hi,H<2\ for the unitary U\$. Following the steps above, we find that 



Vj(e,H u H 2 ) = C/i(0(l + ei))[/i(4ei^cos($)) [l + ie?M 3 ($, H x , -H 3 ) + 0{e\)] 



(A9) 



Assuming e 2 is small, 



cos($) 



cos(</>) cos(0£2) — sin(0) suv{(f)e 2 ) 
cos(0) — sin(0)</>e 2 



(A10) 



and 



Vj{9 7 H u H 2 ) « Ui(9) [1 -ieie 2 4^sin(0) 



^M 3 ($,ffi,-fl- 3 )] • 



(All) 



The key result is that for fixed £2 the infidelity scales as ef when ei is large and e 2 e1 when ei is small. If we 
can improve the U 2 pulses by compensation, we can reduce the error to higher order in e 2 . This is exactly how the 
BB1-WJ sequence (Equation [7]) is constructed, resulting in a fidelity that scales as ef when ei is large and e\e 2 when 
ei is small. 
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